
trtD ■£>*># ® 


2 - 


///3 3 / 
l\ % " 


Bandwidth Efficient Block Codes 


for M-ary PSK Modulation 

i NAS A-CB- 18 1 5 <43 ) BANDWIDTH EFFICIENT BLOCK N88- 12712 

CODES FOB B-ary EBB FGDULATICK {Hawaii 
Uriv.) 48 p CSCE 17B 

Unclas 
G3/32 0111331 


Technical Report 
to 

NASA 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 


Grant Number NAG 5-931 


Shu Lin 

Principal Investigator 
Department of Electrical Engineering 
University of Hawaii at Manoa 
Honolulu, Hawaii 96822 


December 31, 1987 


Bandwidth Efficient Block Codes 
for M-ary PSK Modulation 


Tadao Kasami 
Osaka University 
Toyonaka, Osaka 560, Japan 


Shu Lin 

University of Hawaii 
Honolulu, Hawaii 96822 


ABSTRACT 

In this report, a class of bandwidth efficient block codes for M-ary PSK 
modulation is presented. A soft-decision decoding for this class of codes is 
devised. Some specific short codes for QPSK, 8-PSK and 16-PSK modulations 
are constructed. These codes have good minimum squared Euclidean 
distances and provide 2 to 5.8 dB coding gains over uncoded QPSK modulation 
without (or with little) bandwidth expansion. The complete weight 
distributions of these specific codes are determined. Based on these weight 
distributions, their error probabilities are evaluated. Some of these codes have 
simple trellis structures and hence can be decoded by Viterbi decoding 
algorithm with relatively simple implementation. Moreover the codes are very 
suitable for use as inner codes for various cascaded coding schemes with Reed- 
Solomon codes as outer codes. 


Bandwidth Efficient Block Codes 
for M-ary PSK Modulation 


Recently a great deal of research effort has been expended in bandwidth 
efficient coded modulations for achieving reliable communications on band- 
limited channels [1-27]. This new technique of coded modulation is achieved by 
coding onto an expanded set of channel signals (relative to that needed for 
uncoded transmission). Coded modulation can provide significant coding gain 
over an uncoded system with little or no bandwidth expansion. Most of the 
research works on coded modulation has been focused on trellis coded 
modulations (TCM), i.e., trellis (convolutional) coding with expanded signal 
sets. Not much has been done on block coded modulations. 

In this report, we investigate block coding for M-ary PSK modulation. 
First we present a generalization of coset codes over binary lattices [20] to codes 
over additive groups. From this generalization, a method for constructing 
block coded M-ary PSK codes is devised. Then a soft-decision decoding 
algorithm for these M-ary codes is provided. Some specific QPSK, 8-PSK and 
16 -PSK codes with good minimum squared Euclidean distances and trellis 
structure are constructed. These codes provide 2 to 5.8 dB coding gains over 
uncoded QPSK modulation. Complete weight distributions of these codes are 
derived. Based on these weight distributions, we are able to analyze the error 
performance of these codes for an AWGN channel. Upper bounds on the error 
probabilities are obtained. Since these codes have simple trellis structure, they 
can be decoded with Viterbi decoding. 

In our next report, we will investigate various cascaded coding schemes 
with bandwidth efficient M-ary PSK codes constructed in this report as the 
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inner code. Prel iminar y results show that large coding gains can be achieved 
over the uncoded QPSK modulation. 


2, Code Construction ove r Additive Groups 

Let A be an additive group (finite or infinite) on which a distance 
between two elements s and s', denoted d(s,s'), is defined. The distance 
measure, d(s, s’), satisfies the following conditions: 

1. d(s,s') = d(s-s',0), (21) 

where 0 denotes the zero element in A, and 

2. d(s,s') = 0 if and only if s = s’. (2.2) 

Let B^, B 2 , *** , Bj£ be Jfc nonempty finite subsets of A for which the 
following unique decomposition property holds : for s^ and sj' in Bj, 

s l + s 2 + — + Sj£ = s i + s 2 + *** + s i» (2.3) 

if and only if Sj = s^’ for 1 £ i <> £. Let S be defined as 

S= Bj +B 2 + •■• +B 4 

= { 8j[ + S 2 + *•* + s^ : e Bj with l^i^A}. (2.4) 


Clearly S is a subset of A. 

For a nonempty finite subset B of A, let d[B] denote the minimum 

distance between elements of B. If B has only one element, let d[B] be defined 
as 00 . For 1 < i £ J&, let d^ be defined as follows: 
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dj = d[B| + + *** + 


(2.5) 


For a positive integer n, let A n denote the set of all n-tuples over A. Define the 
distance between two n-tuples s = (s^, 82* *, Sjj) and s’ = (s^ , s£ , *•*, s^) over 

A, denoted d^ n \s, s’), as 

n 

d< n \s,8') = £ d(sj,Sj'). (2.6) 

j=l 

The siim of two n-tuples over A is defined as the component-wise sum of the 
two n-tuples. For 1 ^ i < & , let Cj be a block code of length n over with 

minimum Hamming distance 8p From Cj, , C2> — ,C&, a block code C of 
length n over S is constructed as follows: 


C = { v^+V2 + — + v^:v^6 Cj for 1 ^ i £ A }. 


( 2 . 7 ) 


We will use the following expression for C, 

C = c x + C2 + "* + Cjj. 

Let I X I denote the number of elements in a finite set X. Then 
£ 

I c I = n i Ci I . (2.8) 

i=l 


Lemma 1 provides a lower bound on the minimum distance of the block code C. 
Lemma 1 : The minimum distance of C with respect to 6 S n \ denoted D[C], is 
lower-bounded as follows: 
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( 2 . 9 ) 


D[C] £ Min Sidi 


Proof : For different v and v' in C, let v and v’ be expressed as 


v = vi+V2+ - + V A » Vf € Ci» 

v' = v± + V2 + • + vi , v- € Cj, 


(210) 


where 


v i “ ( ^l* ®i2» » 8 in^ » ®ij e ®i» 

v i = ^ ^l’ s i2» ’ ‘ ’ 8 in^ » ®ij e ®i» 


( 2 . 11 ) 


with 1 < i ^ A and 1 ^ j ^ n. Let h denote the first suffi x such that 


vh* v h- 


( 2 . 12 ) 


Then, since the minimum Hamming distance of is 8^, there exist 8^ 
suffices 1 <J2 < • < js^ - n such that 


s hj p * 8 hj p ’ forl:Sp£5 h- 


( 2 . 13 ) 


Since sy = sjj for 1 <. i < h and 1 S j < n, we have that, for 1 <pi Sjj, 



8 ij + B h + B h+1 + - + B i ] - 


( 2 . 14 ) 
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It follows from (2.1), (2.5) and (2.14) that, for 1 £ p £ 5^ 


% % 

d < I ®ij -X s ii )^ d h 
i=l J p i=l U P 


n 8 8 

Since d^ n \v,v' ) = X d( X 8 ij>.X ty , we have that 
j=l i=l J i=l 


d(n)( v ,v' ) > Sudu S Min Sid; . 

l£i£8 


(2.15) 


(2.16) 


3. Construction of Block Codes for M-arv PSK Modulation 


In this section, we consider code construction for M-ary 

PSK modulation 

where 


M = 2® . 

(31) 

Let 


A = { 0, 1, 2, - , M-l) 

(3.2) 

be the integer group under the modulo-M addition. Define a distance between 

two elements s and s’ in A as follows: 


d(s,s') = 4sin 2 (2" A 7l:(s-s , )) . 

(3.3) 

It is clear that d(s,s') = d(s-s',0) and d(s,s) = 0. For 1 £ i <, let 


Bj = ( 0, 2 i ' 1 }. 

(3.4) 
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t 

Then, B^, B 2 , * , have the unique decomposition property which is related 
to standard binary representation of an nonnegative integer. Note that 

A = B^ + B2 + •” + Bj^, (3.5) 


Then, it follows from (2.5), (3.3) and (3.4) that 

dj = 4sin 2 (2*~^~^rc), for 1 <1 i < &. (3,6) 


Since I B^ I =2 for 1 £ i £ A , a block code over with minimum Hamming 
distance §j can be derived from a binary block code of the same code length 
with the same minimum Hamming distance 5j by substituting 2 1— ^ for 1 in 
each component of a codeword in C^. The code over will be denoted by 
2 i ~ 1 C b . 

Suppose that is a binary linear (n,k^j) code with minimum 
Hamming distance Sj for 1 ^ i £ 9>. Let C denote the sum code, + 2Cfo2 + 
- + 2*-lC b2 . C is linear code over A. C^j is called a binary component code 

of C. It follows from (2,8) that 


ICi 


% 


= 2 i=i kbl 


(3.7) 


Let s and s’ be two n-tuples over the group A. It follows from the definition of 
d (n) (s,s) given by (2.6) and the definition of d(s,s') given by (3.3) that d^ n \s,s') is 
simply a squared Euclidean distance between the two n-tuples s and s’ over A. 
.The minimum squared Euclidean distance (MSED) of code C is then given by 

D[€l = min{S n \vy ) : v, v’e C and v * V ) . (3.8) 
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It follows from Lemma 1, (2.9), (3.6) and (3.8) that 


D[C] 2 Min 4& sm 2 ^' 1 ' 4 *) . 
l£i£4 


(3.9) 


If each component of a code vector v in C is mapped into a point in the 2- 
dimensional 2^-PSK signal set, we obtain a block coded 2^-PSK code. The 
effective rate of this code is given by 


RfC ]=ir| 1 k bi 


(3.10) 


which is the number of information bits transmitted by C per dimension. Let 
C g denote a standard reference code. The asymptotic code gain of C, denoted 

7 (C], over the reference code is given by [8,20] 


, nl R[C]D[C] 
*C] - 101og 10 R[CJD[C a ] 


(3.11) 


If the uncoded QPSK is used as the reference code C s , then R[C S ] = 1, D[C S ] = 2 
and 


tfC] = 101og 10 


R[C]D[C] 

2 


(312) 


The asymptotic coding gain is used as a simple measure of the performance of 
a code. To analyze the performance of a code in details, we need to know the 
complete weight distribution of C. 

Let v = (v| , V 2 , **■ , v n ) be an n-tuple over the group A. The composition 
of v, denoted comp(v), is a M-tuple 
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t 

t = (to.tl* " ^M-l) 

where tj is the number of components vj in v equal to the integer i in A. Let 
W(t) be the number of codewords v in C with comp(v) = t. Let T be the set 

T =* {(tfl, t^, , tj^) : 0 £ tj £ n with 0 £ i < M }. (3.13) 


Then 


(W(t) : t e T } (314) 

is the detail weight distribution of C. Once this weight distribution is known, 
the error performance of C can be analyzed and computed. 

W(t) can be enumerated from the joint weight distribution [28] of the 
binary component codes, 0^1 > ^b 2 > — ,0^^. For a binary A -tuple h = (h^, h 2 , 
*** , hj^) e {0,1}^ and binary vectors v^ = (v^, v^, — , v^ n ) with 1 £ i <, 9, t let 

comp( h ; v^, V 2 , ,v$ ) (315) 

denote the number of j's such that Vjj = hj for 1 £i £ For nonnegative 
integers tQ, t lt ,tM-l» let 


Wj(t 0 , t lf ,t M -i ) 
denote the number of A -tuples, 
( v l» v 2» — 


(316) 
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with Vj € 0 ^ for 1 £ i £ A such that 


comp( h ; v lt v 2 , ••• »v 2 ) = tjj 

for 0 <> h < M, where h is the standard binary representation of integer h. It 
follows from the construction of C that 


W((t 0 , t^, — ,tj^-i)) = Wj(tQ, tp — >t|£_j). (3.17) 


The set, 


* 


{ Wj(t Q , tj, - ,^ 4 ) } 


is the joint weight distribution of the binary component codes, C^i , C\y2’ “* » 
c bi- 

If a maximum likelihood decoding algorithm is used, it is desirable for a 

code to have a simple trellis structure [20]. A trellis diagram of C is a direct 
product of those of binary component codes < 44 , 0 ^ 2 * **• > The num ^ er 

states of a trellis diagram of a binary (njt) code is upper-bounded by [29] 

2 min{k,n-k} ( 3 .i 8 ) 

Some codes may have a trellis diagram with smaller number of states than the 
bound. For instance, the (16,11) Reed-Muller code has a 4-section trellis 
diagram with 8 states [20], The number of states depends on the order of bit 
positions. It can be proved based on Appendix 1 in [ 20 ] that the number of 
states is equal to (3.18) for any n- section trellis diagram of a shortened cyclic 
code. The order of bit positions should be chosen in a clever way. 
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The code construction presented in this section is actually a 
generalization of Sayegh's [26]. 


4, Some Blo ck Codes for QPSK.8-PSK and lfrPSK 

In this section, we construct some block codes for QPSK, 8-PSK and 16- 
PSK modulations. These codes have good minimum squared Euclidean 
distances. Some of these codes have simple trellis structure and hence can be 
decoded by Viterbi decoding algorithm. The codes are constructed in such a 
way that they are suitable for being used as inner codes of various cascaded 
coding schemes with outer codes over GF(2®). For QPSK, the signal set is 

shown in Figure 1. The construction of QPSK codes is based on the additive 
group, A4 = { 0, 1, 2, 3 }, modulo-4 with = { 0, 1 } and B2 = { 0, 2 }. It follows 

from (3,6) that the distances between signal points are : 

d^ = 2, d 2 = 4. (4.1) 

For 8-PSK, the signal set is shown in Figure 2. The symbols for 8-PSK codes 
are from the group, Ag = ( 0, 1, 2, 3, 4, 5, 6, 7 }, modulo-8 addition. = { 0, 1 }, 

B 2 = { 0, 2 } and B3 = { 0, 4 } are chosen to be the symbol sets for the component 
codes. From (3.6), we find the distances between signal points are: 

&l = 0.586, d 2 = 2, d 3 = 4 . (4.2) 

For 16-PSK, the signal set is shown in Figure 3. The 16-PSK codes have 
'symbols from the group, A^g = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 }, 

under modulo-16 addition. The four component codes have symbols from B^ = 
{ 0, 1 } 3 B 2 = { 0, 2 } and B3 = { 0, 4 } and B4 = { 0, 8 } respectively. From (3.6) we 
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find that the distances between signal points are d^ = 0,152, d2 = 0.586, dg = 2 
and d4 = 4. 

Let P n denote the (n,n-l) linear binary code which consists of all the 
even weight vectors. The minimum Hamming distance of P n is 2. Let P^ 

denote the dual code of P n . Then P^ is the (n,l) code which consists of the all- 
zero and all one vectors. The minimum Hamming distance of P^ is n. Let 
^2 m -l denote the even-weight subcode of the binary (2 m -l ,2 m -m-l ) 
Hamming code. Then the minimum distance of is 4. In the following, 

we present a sequence of specific codes for QPSK, 8-PSK and 1 6-PSK. 

Example 1 : Let A - 2, M = 2^ = 4 and n = 5. We choose Cfcj = ?2 x P3 and = 
{0,l} s as the two binary component codes where P2 x P3 is the cartesian 


product of P2 and P3. The binary component code is a (5, 3) code 

with minimum Hamming distance 5^ = 2. The second binary component code 
Cb2 is simply the (5, 5) code. Let 

eg = c bi +2c b 2 - 


(2) 

Then Cq ^ is a QPSK code over the additive group A4 = { 0, 

1, 2 , 3 } with the 

following parameters: 


1 cgj 1 = 28, 

(4.3) 

D[Cq 2 i 1 = 4. 

(4.4) 

r[C q 2 1 ] =I* 

(4.5) 
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(4.6) 


y « 101 og 10 ^ = 2 . 0 (dB). 


This QPSK code maps a 8 -bit message into a sequence of 5 symbols over A 4 . 

Each of these 5 symbols is then mapped into a two-dimensional signal point in 

( 2 ) , . 

the QPSK signal set shown in Figure 1 . The QPSK code Cq ^ can be shown to 
have a trellis diagram with two states. 

FTamplfi 2 : Let A = 2 , M = 2 ^ = 4 and n = 15. Chose C^i = and - ?X 5 be 
the binary component codes. Then the minimum Hamming distances of C^i 
and C ^2 are 8 ^ = 4 and 82 = 2 respectively. Let 


C Q 4 3= C bl+ 2C b2- 


ThenC 


(4) 

Q,3 


is a QPSK code with the following parameters: 


1 eg 1=224, 

(4.7) 

(4) 

D[C^] = 8 , 

(4.8) 

R[C g = 1 ’ 

(4.9) 

y - 101 og i0 jr- = S.CKdB). 

(4.10) 


Since Hi 5 is obtained by truncating the (16,11) Heed-Muller code, the 
component code C^i has a 4-section trellis diagram with 8 states[20]. The 
component code C ^2 bas a trellis of two states. 
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( 4 ) 

The detail weight distribution of the QPSK code Cq g can be evaluated 

easily. For integers i, j and h such that i and j are even, 0£h£i<15 and 0 £ h 
£ j £ 15-i+h, 

W((15-i-j+h, i-h, j-h, h)) = A H>i ^ Jjfh) • (4-11) 

where Ajjj denotes the number of codewords in H-^g with weight i. For other 
composition, t = (tQ, t^, t^, tg), 

W(t) = 0. (4.12) 

In the next 5 examples, we present specific codes for 8-PSK modulation. 

tfoflrapk 3: Let % = 3, M = 2 3 = 8, n « 8. Choose C^i = Pg, 0^2 = ?8 ^ c b3 = 

{0, 1}^ as the 3 binary component codes. The minimum Hamming distances of 
C M , 0^2 and Cj^g are 5^ = 8, 82 = 2 and 63 = 1 respectively. Let 

^8,2 ~ ^bl +2 ^b2 +4 ^b3° 

( 2 ) ( 2 ) 

Then Cg 2 is a 8-PSK code with symbols from Ag » { 0, 1, 2, 3, 4, 5, 6, 7 }. Cg^> 

has the following parameters: 


1 eg I =216 

(4.13) 

D[0g] = 4, 

(4.14) 

Rfcg] = i, 

(4.15) 
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Y = 101og 10 2 = 3(dB). 


(4.16) 


Note that this code provides a 3 dB (asymptotic) coding gain over the uncoded 
QPSK modulation without bandwidth sxs msian- c bl 311,1 C b2 both have trelUs 
diagrams with two states, eg has an 8-section trelUs diagram of 4 states 

as shown in Figure 4 (see Appendix A for construction). 

Hence it can be decoded with Viterbi decoding algorithm. The implementation 
should be rather simple. This block 8-PSK code may be considered to be 
equivalent to Ungerboeck's 4-state trellis code^for 8-PSK modulation which has 
squared minimum free Euclidean Distance d^ = 4 [8]. 

The complete weight distribution of the 8-PSK code C 32 can eva ^ uate ^ 

from the joint weight distribution of its binary component code. For integers 1 , 
j and h such that i are even, 0 < j < i ^ 8 and 0 < h < 8~i, 


W((h, 0, j, 0, 8-i-h, 0, H, 0)) = W((0, h, 0, j, 0, 8-4-h, 0, H» 


'8YiY8— P 

Ai A h )' 


(4.17) 


For other composition t = ( tg, t^, t 2 > tg, 1 4, t5, tg, 1 7 ), 


For its bandwidth efficiency, coding gain and simplicity m 
implementation, eg 18 extremely suitable for use as the inner code for a 

cascaded coding scheme with the NASA standard (255,223) RS code over GF(2 8 ) 
as the outer code. Our preliminary results show that large coding gam can be 
achieved by such a cascaded coding scheme. 
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Example 4 ; Let A =s3, M = 2^ = 8 and n s 23. Consider the 8-PSK code Cg y with 

binary component codes C^, C^2 and where is the (23,12) Golay code 
and Ck 2 =Ck 3 =P 23 . The minimum Hamming distances of C^i , C^2 and C^g 
are 5^ = 7, 52 = 5g = 2 respectively. Then 

C 8?7= C bl +2C b2 +4C b3 

has the following parameters: 


1C® I =256, 

(4.19) 

D[cg] = 4, 

(4.20) 

(2) 28 
Kicy-f. 

(4.21) 

Y = 101og 10 g=3.8(dB). 

(4.22) 


has a 3-section trellis diagram with 2^ states[20j. The complete weight 

( 2 ) 

enumerator of Cg y can be derived from the Hamming weight enumerator of 
the (23,12) Golay code. 


Ftamnlfl fv Let A = 3, M = 2 3 = 8, n = 15. Let C^ , C^2 and be the shortened 
(15,4) Reed- Muller code, P^g and { 0, 1 respectively. Then, = 8, $2 = 2 and 

63 = 1. Let 


C 8,4 " G bl +2C b2 +4C b3- 



(2) 

Then Cg 4 is an 8-PSK code with the following parameters: 


1 eg 1=233, 

(4.23) 

D[cg] = 4, 

(4.24) 

Rrc <®i _ 11 

Klc 8,4 J “ 10 » 

(4.25) 

y=101og 10 y=3.4(dB). 

(4.26) 


Cfci has a 4-section trellis diagram with 8 states[20]. The complete weight 

( 2 ) 

enumerator of Cg 4 can be derived from the Hamming weight enumerator of 
the shortened (15,4) Reed-Muller code. 

Example 6 : Let & = 3, M = 2^ = 8 and n = 4m+3 where 4 £ m £ 7. Let C^, 0^2 
and 0^3 be P^» the shortened (n,n-6) code of H34 and P n , respectively. Then, 

= n, $2 = 4 and 63 = 2. Let 


C 8,m “ G bl +2C b2 +4C b3- 


Then 



is an 8-PSK code with the following parameters: 


IC 8 ( ^l =28“ 
(4) 

«C ] - 8 > 



4m 

4m+3 * 


(4.27) 

(4.28) 

(4.29) 


16 - 



r=ioi°eiorar • 


(4.30) 


(4) 

The complete weight distribution of Cg m is known. For integers i, j and 
h such that i and j are even, 0 < h £ i ^ n and 0 £ h ^ j ^ n-i+h, 


W((n-i-j+h, 0, i-h, 0, j-h, 0, h, 0 )) 


= W(( 0, n-i-j+h, 0, i-h, 0, j-h, 0, h )) 



(4.31) 


where Ajj n j denotes the number of codewords in 0^2 with weight i. For other 
composition t = ( tQ, t^ , t 2 ? tg, t 4 , tg, tg, iq ), 


W(t) = 0. (4.32) 

and Cjjg have trellis-diagrams with two states. Since derived by 

shortening the (32,26) Reed- Muller code, it has a 4-section trellis diagram with 
16 states [20]. 

Example 7 : Let A = 3, M = 2 3 = 8 and n = 29. Let C^x is the linear (29,5) code 

which is obtained from the (32,6) Reed-Muller code of minimum weight 16 by 

first deleting two redundant bits and then truncating one information bit. The 
minimum weight of C^i is at least 14. Let 0^2 b e the linear (29,23) code 

obtained from Hgx by truncating two information bits, and C^g be ^29* ^ben 

§2 ^ 4 and 5g = 2. Let 
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C$ iG bl+2C b2+ 4Cb3- 


Then C 


(4) 

8,7 


is an 8-PSK code with the following parameters: 


i Cg j 1 =2 56 , 

(4.33) 

(4) 

D[Cgj] = 8, 

(4.34) 

«cg]-|. 

(4.35) 

Y=101og 10 ^=5.8(dB). 

(4.36) 


Each of and C^2 has a trellis diagram with 16 states. 

Example 8 : In this example, we construct a code for 16-PSK modulation. The 
code has symbols from the group, A-^g = { 0, 1, 2, ••• , 15 }, modulo-16. The four 

binary component codes, C^, C^3 and (^4, used in the construction are 

p32» the second-order (32,16) Reed-Muller code, P32 and { 0,1 }^, respectively. 

Let 

c 16a0 = C bl +2C b2 +4C b3 +8C b4- 

( 2 ) 

Then C-j^g iq is a 16-PSK code which has the following parameters: 


IClSo' = 2 * 0 ' 


(4.38) 



(4.40) 


* C 16 ioJ = ! 


y = 101og 10 j = 3.9(dB). 


(4.41) 


Since the second-order (32,18) Reed-Muller code has a trellis diagram with 2^ 
states [20], has a trellis diagram with 2^ states, and is invariant to 

phase shifts of all multiplies of id 8. 


5, Encoding and Decoding 

In this section, we consider the encoding and decoding of the M-ary PSK 
codes constructed in Section 3 with M = 2^. Let C denote a M-ary PSK code. 

Encoding 

In encoding, a k-bit message u is divided into 9. submessages, u^ , U 2 , 
• • • , u^ such that the i-th submessage u^ consists of bits and k = k^ + k^2 
+ •*" + kfo^. For 1 ^ i < 9, the i-th submessage Uj is encoded by the binary 
component code Cbj encoder. Let v^ = (v^, Vj 2 ? — , Vj n ) be the codeword for u^. 

Then the codeword for the entire message u is 

v = ( 8J, s 2 , - , Sjj) = V! +2v 2 + - + 2 4_1 V 4 , (5.1) 

where 

sj = V!j + 2v 2 j + + 2 A “" 1 v A j , (5.2) 

for 1 £ j £ n. Note that Sj is a symbol in the additive group, A = {0, 1, — , 2 a -l], 
modulo-2^. The symbols, s^, sg, s Sn are then modulated (mapped into 
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points in the 2-dimensional 2*-PSK signal set) and transmitted. The overall 
encoder is shown in Figure 5. 

Note that the correspondence between the message u and codeword v is 
one-to-one. Let f denote the mapping of u onto v and P* denote the inverse 
mapping of f. Then 


v = f(u) , 

and (5.3) 

i = H(v). 

The mapping f depends on how to divide the message u into £ submessages. 

DmttHifig 

In the following, we present a soft-decision decoding algorithm for the 
M-ary PSK code C. 

For s in { 0, 1, *** , 2^-1 }, let X(s) and Y(s) be defined as 


X(s) = cos(2l~*jts). 

(5.4) 

Y(s) = sin(2 1 - i it3). 

(5.5) 


For s and s’ in { 0, 1, , 2^-1 }, it follows from (3.3), (5.4) and (5.5) that 

d(s,s’) = (X(s>-X(s , ))2 + (Y(s >-Y(s‘)) 2 . (5.6) 

For 1 £ n, let (xj,yj) be the normalized output of the coherent 

demodulator [30] for the j-th symbol of a received vector. The received vector is 
then represented by the following 2n-tuple: 

z = ( y 1? *** > Xjj, y n ). 


~2Q~ 



For the received vector z and a codeword v = (s^ , 82, — , s n ) in C, let I z,v 1 2 
be defined as follows: 

n 

I z,v 1 2 = X (xj-X(sj)) 2 + (y r Y( 8j ))2 (5.7) 


We assume that the channel is an AWGN channel. When symbol s e { 0, 1, 
— , 2*-l } is transmitted, the normalized output (x, y) of a coherent 
demodulator for 2 A -ary PSK is distributed with joint probability density 
function, 


p(x,y) = 


1 -{(x-X(s)) 2 + (y-Y(s)) 2 ]/2c 2 

2X0 2 


(5.8) 


where 



(5.9) 


and p is the SNR per symbol [30, p. 167]. We also assume that every codeword of 
C is transmitted with the same probability. 

Decoding rule : For a received vector z, choose a codeword v in C with 

minimum I z,v I 2 . Then the decoded message u is given by u = f*^(v). 

This decoding rule achieves maximum likelihood decoding for C over an 
AWGN channel. 

If C has a simple trellis structure (the number of states is moderate), the 
decoding of C can be implemented with Viterbi decoding algorithm. 
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Assume that the channel is an AWGN channel and every codeword of C 
is transmitted with the same probability. Let P c be the probability that a 

decoded vector is error-free and Pj c be the probability that a decoded vector is 
erroneous. Since C is linear over { 0, 1, ■■■ , 2^-1 } under addition modulo- 2^ 
addition, we assume that the zero codeword 6 is transmitted without loss of 
generality. For a received vector z, the decoded vector is error-free, if and only 
if 

I z,v I ^ > I z,0 1 ^ , (6.1)* 


for every nonzero codeword v of C. It follows from (2.6), (3,3), (5.4), (5.5) and 
(5.7) that the inequality of (6.1) can be rewritten into the following inequality: 


n n 

(X(sj)-1 Xxj-1 )+Y(sj)yj < £ (X( 3j )-l) 2 +Y(Sj) 2 = d< n >(v,6) . 
j=l j=l 


( 6 . 2 ) 


For an n-tuple v = ( s^, S2 » •** , Sjj) over the group { 0, 1, — , 2^-1 }, 
let Q(v) be the set of vectors, ( x^, y^, — , x Q , y n ), which satisfy the inequality 
of (6.2). Define Q c as follows: 



The probability that I z,v 1 ^ = I z,6 i ^ is zero, and such a case can be neglected. 
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1 

(27ta 2 ) n 




-( I(xj-l) 2 +yj 2 )/2a 2 
e j=i dx^dyj •• dx n dy n , 


(6,4) 


where the integration is taken over Q c , Numerical computation of the integral 
is not feasible unless n is small or Q c has a simple structure. 

The following lemma holds on Q(v). 


Lfiaraa 2' Letv = ( s 1} s 2 , , s^, V = ( s{ y s 2 , ••• , s^)and v" = ( aft, 

s”) be n-tuples over { 0, 1, — , 2^-1 }. Then 

Q(v) n QCv') c Q(v”), 

if the following condition (i) or (ii) holds: 

(i) For each j, one of the following conditions holds. 

(1.1) Sj’ = sj and Sj=0, 

(1.2) sj' = Sj and sj = 0, 

(1.3) sj’ = 2 4 - 1 and ^ = 2 s - 1 + Sj ( mod 2* ). 

(ii) (ii.l) Sj’ = 0 or a' = for 1 <j 2 n, 

(11.2) if 8j = 0 S then sj’ = sj, and 

(11.3) there is an s such that X(s) > 0 and Y(s) £ 0, and if Sj = , then 

either sj = s and 8j' = 2 A ~* - s or sj = 2^ — s and = 2 + s. 

Proof : (i) Suppose that condition (i) holds. Then we have that for 1 < j £ n 

X(sj) - 1 + X(g') - 1 = XCsp - 1, 


Y(sj) + Y(^) = Y(sp , 
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(X(sj) -1)2 + Y(sj) 2 + (X(g') - 1) 2 + Y(^') 2 = (X(sj') - l) 2 + Y(sj')2. 


Hence, inequalities (6.2) for v and v' imply inequality (6.2) for v' ' . 

(ii) Suppose that condition (ii) holds. Then we have that for 1 S j ^ n, 

X(sj) - 1 + X(sXX(^) - 1) = X(sj') - 1, 

Y(sj) + X(s) Y(sj) = Y(sj') , 

(X(8j) -1)2 + Y(sj) 2 + X(s){(X(sj) - 1) 2 + Y(^) 2 ) 

= (X(g") - l) 2 + Y(sj') 2 . 

Hence, inequalities (6,2) for v and v imply inequality (6.2) for v' ’ . ■ ■ 

For a set T of n-tuples over { 0, 1, — , 2^—1 }, a subset S of T is said to be 
T-representative, if 


n 

V6 T 


Q(v) = 


n 

veS 


Q(v) 


(6.5) 


In the examples below, a relatively small subset S of nonzero codewords of C 

can be chosen as a { C — (6} ^representative set. 

For nonzero codeword v of C, let P e (v) denote the probability that a 

received vector z satisfies the following condition: 

I z,v I ^ < I z,6 I ^ , 


that is, 
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n 

2£ (X(sj)-1 Xxj-1 )+Y(Sj)yj S Ivl 2 , 

j=l 


where I v I denotes 



Since the random variable, 


n 

2£ (X(sj)-1 Xxj-1 )+Y(sj)yj , 
j=l 


has a Gaussian distribution with zero mean and variance I v 


p e (v)= J - pi L . e 

,- |2 ^solvl 


x 2 

8o 2 1 v 1 2 dj 



2V5o 



■) 


where 


erfc(x) = 



J e * 2 dt 


and p is the SNR per symbol [30]. 
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( 6 . 6 ) 


I we have 


(6.7) 



For a set Q of n-tuples over { 0, 1, ••• , 2*-l }, let Q denote the 
complementary set of Q. Suppose that S is {C-{6}}-representative. Then it 
follows from (6.3) and (6.5) that 

Hence we have the following upper bound on ^ic* 


Pic = l-P c S I P e (v). (6.8) 

veS 

Let A be the set of real numbers d such that there is a nonzero codeword 

v in C with squared Euclidean distance d from the zero codeword 0. For de A 
and a subset S of C, let A^fS] be the number of codewords of C in S with squared 

Euclidean distance d from the zero codeword. Then, it follows from (6.7) and 
(6.8) that 


p ic £ A d® erfc( 

de A 


(6.9) 


A^[C— (6)] can be computed from the complete weight distribution of C. If we 
can choose a small {C-{6}}-representative set S, A<j[S] may be much smaller 
than A^[C~{6)] except for ’’dominant" d's close to D[C]. 

In the following, we evaluate the error performance of the specific codes 
constructed in Section 4. 

(2) 

Example 1 : Suppose C is the QPSK code Cq ^ . The following subset S of C can 
be easily shown to be (0 — {6}} -representative by using lemma 2. 
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S = { ( S}, 82, * , 85 ) : one component is 2 

and the other components are zero. } 

u { ( s^, S2> 0, 0, 0 ) : si € { 1,3 }, S2 € { 1,3 } } 

u { ( 0, 0, S3, S4, 85 ) : two components of { S3, 84, 85 } are in { 1 ,3 } 

and the remaining one is zero. }. 


It follows from ( 6 . 7 ) and (6.8) that 

p ic 5 Yerfd'Jp). 


( 6 . 10 ) 


Example 2 : Let C be the QPSK code Cq 3 . Note that C^i (= H^g) is a shortened 

code of the second-order (16,11) Reed-Muller code, denoted RM 4 2 . A codeword 

v of a linear code is said to be decomposable if there are two nonzero codeword 
and V 2 in the code such that v = + V 2 and the Hamming weight of v is the 

sum of those of and vg • By using the canonical expressions of codewords of 

the second-order Reed-Muller code [28], it can be shown that any codeword of 
RM 4 2 with weight 8 or 12 is decomposable. It follows from this fact ( Appendix 

B) and Lemma 2 that the following subset S of C is { C— { 6 } ^representative. 
Let I vljj denote the Hamming weight of v. 

S = { v=( 82 , > s l 5 ) : ivl H = 2 and s^ e { 0,2 } for 1 ^i <15 } 

u { y-( si, S 2 » — » S 15 ) : I vl H = 4, 6 or 10, and Sj e { 0,1,3 ) 

for 1 £ i £ 15 and the number of 
occurrences of symbol 3 is even, } 
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u { v=( 8 ^ , 82 , *' * S 15 ) : I v I H = 5, 7 or 11 , and the number of 


occurrences of symbol 2 is one, 
and that of symbol 3 is odd. } . 


It follows from (6,7) and (6.8) that 

P ic } 5 ^erfc(V2p) + 9100erfc(V3p) + 40320erfc(2 Vp) 

+ 43008 erfcCVsp) + 21 5040 erfc( VipT . (6.11 ) 

( 2 ) 

Example 3: Let C be the 8-PSK code Cg 2* Suppose that, for each message (aj , 

a 2> » a 16^ j the bit a l use< ^ as the input to the encoder, the bits ag, ag, 

— , a ^5 are used as the input to the 0^2 encoder and the bits a 2 > a^ • , a^g 
are used as the input to the C^g encoder. The following subset S is {C-{0}}- 

representative, 

S = ( ( s i> S 2 » — , sg ) : one component is 4 and the others are zero. } 

u { ( 32, s 23 — , sg ) : the number of nonzero components is 2, 

and a nonzero component is 2 or 6. } 

u { ( S 2 , S 2 » — , sg ) : 3j is 1 or 7 for 1 £ i £ 8 and the number of 

symbol 7 is even. }. 


it follows from (6.7) and (6.8) that 

P-^ £ 60erfc (Vpj + 64erfc(V2(2-V2)p) . (6.12) 
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On the other hand, the probability Pq^ c that there occurs at least one bit error 
when 16 bits are transmitted by uncoded QPSK is given by [30], 


p Q,ic = 1 - [ 1 - | erfc < VlT )] 16 - (6.13) 

( 1 ) 

In Table 1 and Figure 6, we show the upper bounds on P- given by 

(6.12) and Pq j c given by (6.13) for various SNR per information bit, p/2. We 

(2) 

see that the code Cg ^ achieves a 3dB real coding gain over the uncoded 
QPSK without bandwidth expansion at 10 block error rate. In Table 2 and 
Figure 7, we show the upper bounds on the decoded bit error probability 
for Example 3. 

Let T be a subset of { 0, 1, — , 2 a -l }. For v e C, define 


Occ( T,v ) = (the number of occurrences of symbol s in v ) . 
seT 


(6.14) 


Example 4 : Let C be the S-PSK code Cg j. Then the following subset S of C — (6) 
is easily shown to be {C-{0))-representative by using lemma 2. 


S = { v=( S| , S2, — , S23 ) : I v ! H = 2 and Sj e { 0, 2, 6 } for 1 < i < 23 } 


u { v=( s^ , 82, • * , S23 ) : Occ( { j },v ) £ 1 for each je { 2, 3, 4, 5, 6 } 

and Occ( (2, 3, 4, 5, 6} ,v) < 2. } . 


Example 5 : Let C be the S-PSK code Cg 4. The following subset S of C can be 
easily shown to be {C — {0}) -representative by using lemma 2. 
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S = { v=( 8^, &2» **’ » ®15 ) : ! v I h = 1 and ®i € { 0, 4 } for 1 £ i £ 15 } 

u { v=( sj, 82, — , Si5 ) : ivl H = 2andsj€ { 0, 2, 6 } for 1 £ 15 } 

u { v=( s^, 32, ••* , Big ) : Occ( U ),v ) ^ 1 for each je { 2, 3, 4, 5, 6 } 

and Occ({2,3,4,5,6},v) ^ 2. } . 

(2) 

Example 8 : Let C be the 8-PSK code Cg m . The following subset S of C can be 
easily shown to be {C — {0}} -representative by using lemma 2. 

S = { v=( s x , S2, — , s n ) : I v ! H = 2 and Sj € { 0, 4 } for 1 < i ^ n } 

u { v-( si, S2 , , Sjj ) : Sj € { 1, 3, 5, 7 } for 1 ^ i < n and 

Occ( { 5,7 }, v ) is even. } 

u { v-( si , S2> ••• , ^ ) : lvl H = 4, 6 or 10< 1 v! H <n , 

Sj € { 0, 2, 6 } for 1 < i ^ n and Occ< {6}, v ) is even. } 

u{ v=( si, S2, , s n ) : I v I H = 5, 7 or 11 < I v! H <n , 

Sj e { 0, 2, 4, 6 } for 1 £ i ^ n , 

Occ{ {4}, v ) = 1 and Occ( {6}, v ) is odd. } . 


The error performance and coding gains of the codes given in the above 
examples are now being computed and will be included in our next report. 
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7. Conclusion 

In this report, we have presented a class of bandwidth efficient block 
codes for M-ary PSK modulation. A soft-decision decoding for this class of 
codes is devised. Some specific codes with good squared minimum Euclidean 
distance are constructed. The complete weight distributions of these codes are 
determined. Their error probabilities are evaluated. Some specific codes have 
simple trellis structures and hence can be decoded by Viterbi algorithm. Some 
of these codes are suitable for use as the inner codes of cascaded coding scheme 
with Reed-Solomon codes over GF(2^) as outer codes. Our preliminary results 
show that such cascaded coding schemes provide extremely high reliability 
and large coding gains. 

In our next report, we will present the coding gains of the codes 
presented in this report over the uncoded QPSK modulation. 
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Appendix A 


Trellis Diagram for the 8-PSK Code Cg % 

(2) 

The 8-PSK code Cg 2 consists of three binary component codes , C^2 
and Cgg which are Pg, Pg and { 0, 1 }** respectively. Let u be a 16-bit message to 

be encoded. Divide u into three submessages , U2 and ug where U4 consists 
of only one bit, U2 consists of seven bits and ug consists of eight bits. Then u^, 
U2 and ug are encoded based on C^i, Cg2 a^d C^g respectively. Let 

a = ( a^, a2, ag, a 4, ag, ag, 37, ag ) 

b = ( h[ , b2, bg, b4, bg, bg, b7, bg ) 

c = ( cj, C2, eg, C4, eg, eg, 07, eg ) 

be their corresponding binary codewords. Note that a is either the all-zero 
vector or the all-one vector. The codeword b has even weight. 

For 1 ^ £ <* 8, the input to the signal selector of the overall encoder- 
modulator at the &-th time unit is the triplet ( a^, bj£, c$ ). If a# = 0, then ( b$, 
C£ ) selects a signal point from the QPSK signal set shown in Figure 2b. If a$ 
= 1 , then ( b$ , c$ ) selects a point from the QPSK signal set shown in Figure 2c. 

Hence the system switches between two QPSK signal sets. To construct the 
trellis diagram for Cg we need to define the states of the overall encoder- 

modulator. Let ( b^, bg, , b$ ) denote the &-bit prefix of codeword b. Let W( 
bi» b2, — , b& ) denote the Hamming weight of ( b^, b2» — » b$ ), At the £-th 
time unit, the state of the encoder-modulator depends on the bit a^ and the 
number W( b^, b2, * , b^ ). Define the following states: 
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(1) A e represents the states that a* = 0 and W( b^ , b 2 » , b$ ) is even; 

(2) Aq represents the states that a$ = 0 and W( b^, b 2 , , b$ ) is odd; 

(3) B e represents the states that a$ = 1 and W( b}, b 2 > - , b$ ) is even; and 

(4) B 0 represents the states that a$ = 1 and W( bj, b 2 » — , b$ ) is odd. 

Assume that the encoder-modulator starts from the state A 0 at the time % = 0. 

(2) 

Then the trellis diagram for Cg ^ can be constructed easily as shown in Figure 

4. There are two parallel branches (or transitions) between the transition of 
two states; they correspond to c& = 0 and c^ = 1 respectively. 

The encoding of message u is equivalent to tracing a path in the trellis 
diagram. The codeword corresponding to u is a sequence of QPSK signal 
points either from the set shown in Figure 2b or from the set shown in Figure 
2c. 
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Appendix B 


Consider the second-order (16,11) Reed-MuUer code RM4 2 • We use a 
boolean function b(v) for expressing a codeword v in RM42 . 

By using an affine transformation* a codeword v in RM4 2 with weight 8 

or 12 can be represented as one of the following forms [28, p.438]: 

1) If 1 v I jj = 8, then b(v) = X3X2+X3 or b(v) = X3. 

2) If I v I jj = 12, then b(v) = x^X2+l . 

Let y be x^, X2 or xj + X2 . Then the degree of yb(v) is at most 2. Let and V2 

denote the codewords represented by yb(v) and (y+l)b(v), respectively. Then, v = 
V| + V2» v x * 6, V2 * 6, and I v I H = I I h + • v 2 * H • is » v is decomposable. 

Let u be a codeword of C^2 ( = P15 ) such that 
I v+2u I H = I v I jj • 

Let I u I jj y =a denote the number of nonzero components of u in the bit- 
positions for which y=a. If I u I H x _q and I u I H x _q are odd, then I u I H>x _x 

and I ix 1 Hj x 2 =1 are odd - Then ,u, H,x 1 +x 2 =0 and lu, H,x 1 +x 2 =l are even ' 
Therefore we can choose one of x^, X2 and x^ + x 2 as y in such a way that 
I u I H,y=0 and * u ^ H,y=l are even * 

S S 
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Figure 1 


QPSK signal set and squared Euclidean 
distances between signal points. 





Figure 3 


16-PSK signal set and squared Euclidean 
distances between signal points. 



( 2 ) 

Figure 4 4-state trellis diagram for the 8-PSK code Cg 2 







Table 1 . Upper bounds on the probability of incorrect 
decoding for a block with 16 information bits 
for Example 3 and uncoded QPSK 









































Table 2 Upper bounds on the decoded bit 
error probability for Example 3 


SNR per 
symbol 

(dB) 

SNR per 
information 
bit (dB) 

mm 

7.0 

3.98 

4.31E-02 

7.5 

4.48 

1.99E-02 

8.0 

4.99 

8.42E-03 

8.o 

5.45 

3.66E-03 

9.0 

5.97 

1.29E-03 

9.5 

6.49 

4.14E-04 

10.0 

6.99 

1.21E-04 

10.5 

7.49 

3.21 E-05 

11.0 

7.96 

7.78E-06 

11.5 

8.49 

1.39E-06 

12.0 

8.98 

2.52E-07 




































